Driving non-Gaussian to Gaussian states with linear optics 
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We introduce a protocol that maps finite-dimensional pure input states onto approximately Gaussian states in 
an iterative procedure. This protocol can be used to distill highly entangled bi-partite Gaussian states from a 
supply of weakly entangled pure Gaussian states. The entire procedure requires only the use of passive optical 
elements and photon detectors that solely distinguish between the presence and absence of photons. 
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I. INTRODUCTION 

Gaussian entangled states may be prepared quite simply in 
optical systems: one only has to mix a pure squeezed state 
with a vacuum state at a beam splitter, both of which are spe- 
cial instances of Gaussian states in systems with canonical 
coordinates Lllj]. The beam splitter acts as a Gaussian uni- 
tary operation which modifies the quantum state, but does not 
alter the Gaussian character of the state. The resulting pure 
state is a two-mode squeezed state. This state may be used as 
the resource for protocols in quantum information processing. 
In fact, teleportation |3], dense coding |4] and cryptographic 
schemes |3|] on the basis of such two-mode squeezed states 
have been either studied theoretically or already experimen- 
tally realised. For the theory of quantum information process- 
ing in systems with canonical degrees Gaussian states play a 
role closely analogous to that of entangled states of qubits, for 
which most of the theory of quantum information processing 
has been developed. 

However, there are significant limits to what accuracy 
highly entangled two-mode squeezed states may be prepared 
and distributed over large distances. Firstly, the degree of 
single-mode squeezing that can be achieved limits the degree 
of two-mode squeezing of the resulting state. Secondly, de- 
coherence is unavoidable in the transmission of states through 
fibres, and the original highly entangled state will deteriorate 
into a very weakly entangled mixed Gaussian state OJ. For 
finite-dimensional systems, it has been one of the key obser- 
vations that in fact, from weakly entangled states one can ob- 
tain highly entangled states by means of local quantum oper- 
ations supported by classical communication |7] at the price 
of starting from a large number of weakly entangled systems 
but ending with a smaller number of more strongly entan- 
gled systems. The term entanglement distillation has been 
coined for such procedures. Importantly, such methods func- 
tion also as the basis for security proofs of quantum crypto- 
graphic schemes |9]. 

It was generally expected that an analogous procedure 
should exists for the distillation of Gaussian states by means 
of local Gaussian operations and classical communication 



only. Surprisingly however, it was recently proven that this 
is not the case Hot Hill . For example, no matter how the 
local Gaussian quantum operations are chosen, one cannot 
map a large number of weakly entangled two-mode squeezed 
states onto a single highly entangled Gaussian state. Gaussian 
quantum operations [ 5Tq1 fllL fl2ll correspond in optical systems 
to the application of optical elements such as beam splitters, 
phase shifts and x -squeezers, together with homodyne de- 
tection. All these operations are, to some degree of accuracy, 
experimentally accessible. With non-Gaussian quantum oper- 
ations, in turn, one can distill finite-dimensional states out of a 
supply of Gaussian states but the resulting states are not 
Gaussian, and the experimental implementation of the known 
protocols constitutes a significant challenge. 
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FIG. 1: A single step of the protocol. Two pairs of entangled two- 
mode states are mixed locally at 50:50 beam splitters and absence 
or presence of photons is detected in one of the output arms on both 
sides. 



One may be tempted to think that this observation renders 
all attempts to increase the degree of entanglement in Gaus- 
sian states impossible. In this article, however, we discuss 
the possibility of obtaining a Gaussian state with arbitrarily 
high fidelity from a supply of non-Gaussian states employ- 
ing only Gaussian operations, namely linear optical elements 
and projections onto the vacuum. We describe a protocol that 
prepares approximate Gaussian states from a supply of non- 
Gaussian states, which shall be called 'Gaussification' from 
now on. The non-Gaussian states that we use could in par- 
ticular be obtained from weakly two-mode squeezed vacua, 
by the application of a beam splitter and a photon detector. 
Together with this step, the proposed procedure offers a com- 
plete distillation procedure of Gaussian states to (almost ex- 
act) Gaussian states, but via non-Gaussian territory. It is im- 
portant to note that the protocol introduced below is by no 
means restricted to a bi-partite setting. The bi-partite case is 
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only practically the most important one, as it allows in effect 
for distillation of Gaussian states with non-Gaussian opera- 
tions. But this method can, in particular, also be used in a 
mono-partite setting to approximately obtain a Gaussian state 
from a supply of unknown non-Gaussian states. 

The paper is organised as follows: First, we will describe 
the protocol, that generates Gaussian states from a supply of 
non-Gaussian states. This protocol requires only passive op- 
tical elements and photon detectors that can distinguish be- 
tween the absence or presence of photons but that do not de- 
termine their exact number. We then proceed by discussing 
the effect of the protocol in more detail. We will discuss the 
special case of pure states in Schmidt form as well as gen- 
eral pure states. The fixed points of the iteration map will be 
identified as pure Gaussian states and a proof of convergence 
will be given. Finally, we will discuss the feasible preparation 
of finite-dimensional states from a supply of pure Gaussian 
states. 



II. THE PROTOCOL 

The protocol is very simple indeed. We start with a sup- 
ply of identically prepared bi-partite non-Gaussian states. The 
overall protocol then amounts to an iteration of the following 
basic steps: 

1 . The states will be mixed pairwise locally at 50:50 beam 
splitters (see Fig.0. 

2. On one of the outputs of each beam splitter a photon de- 
tector distinguishes between the absence and presence 
of photons. It should be noted that we do not require 
photon counters that can discriminate between different 
photon numbers. 

3. In case of absence of photons at both detectors for a par- 
ticular pair one keeps the remaining modes as an input 
for the next iteration, otherwise the state is discarded. 

This is one iteration of the protocol which we will con- 
tinue until we finally end up with a small number of states 
that closely resemble Gaussian states. This is clearly a proba- 
bilistic protocol. However, the success probability, as we will 
see later, can be quite high. It should also be noted that the 
operations in a successful run are indeed Gaussian operations, 
namely the use of linear optical elements and vacuum projec- 
tions. Each of these steps can be realised with present-day 
technology. 



III. EXAMPLES OF THE PROTOCOL 

A. Pure states in Schmidt form 

In order to demonstrate the general mechanism, we start 
by discussing a particularly simple case, namely pure states 
in Schmidt form. We do not require any prior knowledge of 



the actual un-normalised state vectors except that they can be 
expressed in the following form 



(o)\ 



E a $>> n > 

71 = 



(1) 



where {a4 n}„ = „ 



with anli 



> are proportional to the real 
Schmidt coefficients of the state vector, and {\n) : n 6 N} 

denotes the Fock basis. We only assume a °Q > and it is 
then convenient to consider un-normalized states for which 
we set a °Q = 1. The un-normalized states arising in later 

steps i = 1, 2, . . . are characterised by coefficients { 
These coefficients then become identical to the Schmidt co- 
efficients only after appropriate normalisation. Starting from 
two identical copies of state vectors that have been obtained 
in the ith step of the protocol , i.e. 



|V> W )|V> ( ^ 



(2) 



one obtains after application of the 50:50 beam splitters the 
state vector (L/i2 O Ui2)\t^}\^ (l) )- Here, the beam splitter 



(3) 



is described by (see, e.g., 1 15]) 



where U\2 acts on the amplitude operators of the field modes 
as 



U- 
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-R* T* ) \ a 2 
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where we set T = R = 1/^/2. The resulting un-normalised 
state vector, conditional on vacuum outcomes in both detec- 
tors, is given by 



(0,0|([/ 12 ®{7 12 )|V>| 



E 

71 = 
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where 



a 



(i+l) 



E 

r=0 



(0 



(5) 



(6) 



for n = 0, 1, . . ..The probability of vacuum outcomes being 
detected in both modes is (V( l+1 )|V (l+1) )/|(V W IV (!;) )| 2 - The 
protocol is a Gaussian quantum operation, in the sense that it 
is a completely positive map that maps all Gaussian states onto 
Gaussian states. The interesting feature is that by repeated 
application it also maps non-Gaussian states arbitrarily close 
to Gaussian states, as will be demonstrated below. 

In effect, in each iteration one maps one sequence of co- 
efficients aM = {an^n}n° = o onto another sequence = 
{a?i,7i 1 ' > } r ^Lo' defining the map $ via 



= : $(a 



(7) 
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In the following we use the notation $W = $ and $( 4+1 ) = 
$o$W fori = 0,1,.... The main observation is that in fact, 

provided a^l < cxqq, the sequence of coefficients {a^}^l 1 
converges to a distribution corresponding to a Gaussian state, 
in this special case a two-mode squeezed vacuum. 

In other words, although the initial state was not Gaussian, 
but say, a state corresponding to a finite-dimensional state vec- 
tor of the form 



|^°)> = |0,0> 



(8) 



where ct^°] 6 [0, 1), after a number of steps the resulting state 
is Gaussian to a high degree of accuracy. We will first show 
that this convergence is a general feature of this protocol, and 
we will then discuss the consequences. We start by demon- 
strating that those distributions associated with pure Gaussian 
states are fixed points of the map <£>. 

Proposition 1. - The distributions a — {a n .„j^o of the 

form 



= X 7 ' 



(9) 



A > 0, corresponding to two-mode squeezed states, are the 
only fixed points of the map $. 

Proof. This can be immediately derived from the definition 
of Let us assume that 



= $(a) 



(10) 



holds. It can be verified by substitution that a n . n = A™ is a 
solution of this equation. The uniqueness of this solution can 
be verified by observing that Eq. JlOi also implies ao,o = ot-o o> 
that is, ao,o = 1- Then a\ : \ is a free parameter and once set 
(i.e. as ai i = A) the remaining coefficients are uniquely 



determined. 

These coefficients, for A S [0, 1), in turn correspond ex 
actly to two-mode pure Gaussian states. If A lies outside this 
range, the state is not normalizable. The next Proposition 
states that those distributions associated with Gaussian states 
are not only fixed points of the map but provided ^ 0, 
each sequence of coefficients converges to such a fixed point. 

Proposition 2. - Let = {a4°n}^L w * tn a oo = 1 an d 
< < 1. Then 



lim a (i) = a (oo) 



(ID 



for all n = 0, 1, where a^ 00 ) is a distribution of the type of 
Proposition 1. 

Proof. As before, let us set aW := $M(a( 0) ) for i — 
1,2, The first step is to see that 



a 



(i+l) 



Q' 



1,1 _ „(0) 



a, 



(i+l) (i) 



(12) 



0,0 



o.o 



for alli = 0, 1, Let us first assume that a'f l > 0. Then, as 

can be seen from the definition of $, 

^ 1) ^\ = l(^ + « i >$)4f- d3) 



Hence, as \ = a± \ > for all i = 0, 1, . 
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i.i 



* (0) 
n,i- 



(14) 



Now let us assume that already a n _ x „_i 
0,1, ...and 



> for all i 



lim (i) 



(*) 



(15) 



for some n = l, 2, 

(i+l) 

a n+l,n+l 



(i+l 



. Then, from 

lZjr=O t,r , ra n-r+l,ii-r+l\ r J 
Z^r=0 " r , rLt n— r.n — ■< 



(16) 



n—r,n—r \rJ 



it follows after a few steps that a 
and 



« 



> Oforalli = 0,1,. 



lim 

i — >oo 



lim 



(i+l) 

n+l,n+l 



(i+l) 



1 



which means that 



za n+l, 



(i) 



lim 

i — too 



(i+l) 



(2" +1 - 2)afl 



/ 0) 



(17) 



(18) 



Hence, by induction we find that the ratios of a n 



(i) 



O-n^n converge to the ratio of < a^j < 1 and a U Q = 1 as i — > 
oo. This means that the coefficients correspond to a Gaussian 
□ state as specified in Proposition 1 . In case that afl = an 



(0) 



+ l,n+l 



and 



„(°) 



(i) 

analogous argument can be applied in order to arrive at a = 
1 for alH = 0, 1, ... and 



lim ai'l = 



(19) 



for all n = l, 2, j-j 

This shows formally that the (pointwise) convergence to an 
effectively Gaussian state is generic 0. Putting aside the 



restriction that a 
more detail. 



(0) 
0,0 



1, three cases shall be discussed in 



1. If a[) Q > and < Qqq, then the states converge to 
a Gaussian state. 

2. A special instance is when a °Q > 0, but = 0. Then 
the states converge to a Gaussian state, but to the prod- 
uct of two vacua. 

3. If a °g < a.f\, then the sequence does not converge to 
a sequence of coefficients corresponding to a Gaussian 
state. In particular, this is always the case when 



a n=o. 



(20) 



This follows immediately from Eq. (|6jl as a = for 
all i. 
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In practice, one can actually expect a state that is very close 
to a Gaussian state already after a very small number of steps, 
say, three or four steps. As has already been mentioned, the 
whole scheme is probabilistic. That is, the success probabil- 
ity of actually obtaining the desired state is always less than 
one. In Fig.[2]we show the total probability of success, PsuLess> 
and in Fig.[5]the corresponding fidelity i.e. the overlap 
with the Gaussian state to which the protocol converges, after 
i = 1,2,3 iteration steps. Here, we started with coefficients 



M - 



1 and a[ °\ — X. We see that for a large range of val- 
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FIG. 2: Success probability process after i = 1 (dotted line), i = 2 
(dashed line) and i = 3 (solid line) iteration steps, where the initial 
states werecx |0,0) + A 1, 1). 



where otm]n € (D for all n, to. If the procedure described in 
Sec. |Il]is carried out, using 50:50 beam splitters with appro- 
priate phases, such that T = R=1/ \/2, then, for a large class 
of input states, after repeated iterations of the protocol, a state 
closely approximating a Gaussian state will be obtained. If 
the identical retained states after i iterations of the procedure 
are labeled 



(22) 



we can describe each iteration in terms of the following recur- 
rence relation, 



^ ) = 2-^ i EE(- 1 ) (m+ 



n) — (r+s) 



r=0 s=0 



^r.s^m—r.n — s 



1/2 



(23) 



where again 



(24) 



with a« = {a^ ) m}~ m=0 for i = °> h ■■■■ We will in the 
following write 



n (oo) := lim 



(25) 
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whenever this limit exists. The fixed points of $, charac- 
terised by Om,n € (D, correspond to states which are un- 
changed by one or more iterations of the procedure, and sat- 
isfy $(a(°°)) = a(°°), thus 



(m+n) — fr+s) 



x a (oo) a (oo) 



: : 



for all n, to. We immediately see that 



-,1/2 



(26) 



FIG. 3: Fidelity of the approximately Gaussian state after i = l 
(dotted line), i — 2 (dashed line) and i = 3 (solid line) iterations 
where the initial states were oc |0, 0) + A 1 1 , 1). 

ues for A the fidelity is just below unity, and for A = 0.5 the 
probability of success is still above 0.5. 



B. General pure states 

Suppose now we have a supply of pure states with state 
vectors of the general form 



(o)\ 



(21) 



m.7i— 



*0,0 " 



-H7) 



and thus a 



(oo) 
0,0 



1. (The other possibility, a 



(oo) 
0,0 



(27) 



leads to 



the trivial solution a 



(oo) 



for all to, n.) We also find that 



the coefficients and are the only free param- 

eters. When these values are specified, all other coefficients 
are determined. The general solution of Eq. j26l > is 



(oo) 
*2m,2n+l 



- a (oo) - 

u 2m+l,2n ~ u : 



(28) 



7i 2 I (7i/2)' 



E 



0<s<m;s<n _ 



(72/2) r 



(2s)! (to-s)! (n-s)! 



(29) 



5 



(oo) 
*2m+l,2n+l 



£ 



V(2m+l)!V(2n + l)! 

(71/2)"^ (72/2)' 



0<s<?7i;s<n 



712 



(2s + l)! (m-s)! (n-s)! 



(30) 



where the coefficients 71 , 72 and 712 are usefully expressed as 
elements of the symmetric 2x2 matrix 



7i 7i2 
712 72 



(31) 



and are determined uniquely by the free parameters 
a 2^ > a< 02 an d c 4°i' ) - A specific form for this correspon- 
dence is given in Proposition 4. The coefficients ai^ deter- 
mine an un-normalized state vector |^(r)). In the Fock state 
representation this state vector is given by 



|^(r)) = Q(r)|o,o), 



(32) 



where the operator Q(T) is expressed in terms of T and the 



vector d 1 = [a\, a\) T as 



Q(T) = cxp 



-(dt)^r(dt) 



(33) 



The state vectors |^>(r)) are not normalized, and the require- 
ment that they be normalizable, i.e. (ip(T)\ip(T)} is finite, 
places a restriction on I\ The following proposition takes its 
most concise form when we use the spectral norm which is 
defined as 1 16] 



I l-X'l |oo — vAmax 

where A max is the largest eigenvalue of XX'. 



(34) 



Proposition 3. - If and only if \\T Woo < 1, then \ip{T)) := 
Q(r)|0, 0) is normalizable and represents a pure Gaussian 
state. 



Proof: The matrix T in Eq. (13 1 1 is a complex symmetric 
2 x 2-matrix. Following Takagi's Lemma llfll . there exists a 
unitary matrix U such that 

U T TU = :A, (35) 

where A is a diagonal matrix the entries of which are the 
eigenvalues of vTTT. With b:— Ua we have 



|V(r)) = cxp 



i(6t)*A(6t) 



|0,0) 



(36) 



Because the b\ and 62 commute, this is a tensor product of 
two single-mode Gaussian states. It is now straightforward 
to show that the single mode state vectors are normalizable if 
and only if both diagonal elements of A are smaller than one. 
Then, each of the modes is in a single-mode squeezed state 
lH/ll . The transformation a 1 — ► Ua represents a beam-splitter 
transformation mapping the original modes a onto the modes 
b, i.e. it is a passive transformation. Hence, the resulting state 



vector Eq. d32i is also normalizable. 



□ 



In fact, as can be shown, the state vector Q(T) | 0, 0) is, apart 
from normalisation, equal to the state vector of the two-mode 
squeezed vacuum state S(Z)\0, 0), where 



S(Z) = cxp 



i(dtfZ(dt)-i(dfZt(a) 



(37) 



S{Z) is a generalized two-mode squeezing operator fvh . 
where Z = arctanh(rr)e ier with the polar decomposition 



T = r v e l ° r 



(39) 



Proposition 4. - Suppose we are given a supply of iden- 
tical two-mode pure states with state vectors \ip^) = 
J2 m ,n am,n\m, n), and let 



V2/3 2 ,o ~ 01 A.i-jSi.oAm 
- Pi,oPo,i V2/3 ,2 - ffi A 

where f3 m ,n-= otm, ] n/ a ofi- ^ll r lloo<l then 

i — ^00 ' 

for all n,m = 0, 1, where 

a m , n := (m,n|Q(r)|0,0). 



(40) 



(41) 



(42) 



Proof. To make the proof simpler, we shall use q;q q = 1 
as above. This is merely a change of normalization and does 
not alter the general validity of the argument. Before proving 
the convergence of all coefficients otm,n under 4> to the fixed 
point a^l, as i — ► 00, let us first show that a certain subset 
of coefficients actually reach their final value after a single 
iteration of $. 

The coefficients a 2 m + i 2n a 2m 2«+i reacn zero, their 
fixed point, after a single iteration corresponding to i = 1, for 
all to, n. To see this, note that in the following equation, 



(m+n) — (r+s) 



x a^a {0) 

r,s m—r,n — s 



TO \ / n 

r As 



1/2 



(43) 



renaming the summation indices (r, s) 1— > (to— r, n—s), yields 
an identical sum except for an overall factor of (— l) m+ ™. 
Consequently, for odd values of m + n the whole sum must 

vanish and coefficients of the form a 2 m+i 2n anc ^ a< 2m 2n+i 
vanish after a single iteration step. As a consequence of this, 

the coefficients ct[\, and \ also do not change after 



one iteration. For example, 

(i+l) 



a 



(0 
1,1 



On a 



0,1"!, 



(44) 
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for all i = 0, 1, .... Similarly, a^o anc ^ a< ol a ^ so assume their 
respective fixed points after the first iteration, and thus the 
matrix T is determined to be as in Eq. ( I40l i. 

Now let us show that all coefficients a.m,n do indeed con- 
verge to their respective fixed points a^L as i — > oo. The 
recurrence relations in Eq. d23t can be re-written as 



»(<+!) 



2 — 5" 



2 <«+ EE (-i) (m+n) - (r 



x 




lV2 n 



(45) 



(i) 

Let us assume that all coefficients a^, where r <m, s <n 
but r + s < m + n, do converge to the fixed points a£°^ as 
i — > oo. Then 

Urn c£t 1} = 2 {1 -^ lim afi? 



m . n 



m n 



(r,a)#(0,0)#(m,n) 

:)(:: 



(46) 



x a(°°)a (oo) 



1/2 



Now let us use the substitution <J^ n := a£,n- and we 
obtain, using Eq. (|26li, 

lim ^t 1 ) = 2( 1 -^) lim . 



We see that 6m,n converges to zero as long as 



2 (i- 



< 1, 



(47) 



(48) 



which is the case whenever m + n > 2. However, since 

U) 

we have already shown that all coefficients a m \n, where 

/ o • (t) (i) (i) (i) (t) . (i) ' 
m + n < 2, i.e. ctg q> ' x , g, a\ \, a 2 and , converge 

to a final value after a single iteration, the convergence of all 

other coefficients follows by induction. Note that whenever 

1 1 r | |ao > 1, although the coefficients individually converge 

to their respective fixed points, the state as a whole does not, 

since Q(T) |0, 0) is not a normalizable state vector. rj 



IV. GENERATION OF THE INITIAL STATES FROM 
GAUSSIAN STATES 

So far we did not specify where the supply of initial states 
should come from. In fact, one could use two (weakly) en- 
tangled Gaussian states and feed them into one of the iter- 
ation components shown in Fig. ^ Then, instead of retain- 
ing the state in the case of measuring the vacuum, we now 



retain the state whenever any nonzero photon number is ob- 
tained. Again, only detectors that distinguish between absence 
or presence of photons are needed. Let us start with a sup- 
ply of two-mode squeezed vacuum states the state vectors of 
which can be written in Schmidt basis as 



IV^Hv/T^E^k")' 



(49) 



with q £ [0, 1). 

In general, it will be easier to generate two-mode squeezed 
states with low values of q in an experiment, and using the 
following simple protocol one can use a supply of such states 
to generate a supply of non-Gaussian states which, when used 
as the input of the procedure described in Sec.|n] lead to the 
generation of two-mode squeezed states with much higher q. 

Let us feed two copies of the state of the form as in Eq. ( I49> 
with q <C 1 into the device schematically depicted in Fig. ^ 
and retain those outcomes that correspond to a 'click' in both 
detectors. It does not matter how many photons have been 
measured, and we do not assume that a different classical sig- 
nal is associated with different photon numbers. The projec- 
tion operator fl8il describing this process is 



p = (i-|o)<o|)®(i-|o)<o|). 



(50) 



Although the vacuum projection (as well as the identity op- 
eration) are Gaussian, the difference of them is not, and, 
indeed, we find that when the states used in the protocol 
have sufficiently small q, then this projection approximates 
|1)(1| eg) |1)(1| with high accuracy. Thus, we are not in the 
situation as in Refs. 1 10, Acting with ( 15 0> on two copies 
of the state i49t after rotating them at the beam splitters gives 
the non-Gaussian state with un-normalised state vector 



\9{q;T A ,R A ;T B ,R B )) := 

P \U 12 {T A , R A ) ® il 12 (T B , R B )] |V 9 ) 02 ■ 



(51) 



where again 

U 12 (T,R) = T fll e- R ' & 2 &1 e m2& iT~ A2 (52) 

wdT A ,T B ,R A ,R B £ € with 

|T A | 2 + |J2 A | 2 = |T B | 2 + |JiB| 2 = l. (53) 

For simplicity of notation, let 

ui{q;T A ,R A ;T B ,R B ):= (54) 
tr M [|tt(g; T A , R A - T B , Rb)) (*(g; T A , R A ; T B , R B )\] 
trp(<7; T A , R A ; T B , R B )) (q; T A , R A ; T B , R B )\] 

be the normalised state after application of the beam splitters 
and the two projections, where \im is the partial trace over the 
measured modes. The most appropriate choice for the reflec- 
tivities and transmittivities clearly depends on the value of q 
and on the figure of merit of how one quantifies the quality 
of the output state. However, when q £ [0, 1) is very small, 
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the output state can be made arbitrarily close to a maximally 
entangled state 

p+ = ±= [|0, 0) + e-**|l, 1)] [(0, 0| + e*(l, 1|] (55) 

in 2 x 2 dimensions. Where the phase e 1 ^ depends on the 
phases of T and R in the beam splitter chosen. More precisely, 



lim \\u>(q;t{q),r{q);0, 1) - p + L = 0, 



where 



l*(9)l 



7 2U/2 



\r(q)\ := [l - \t{q)\ 



21 2 



(56) 

(57) 
(58) 



and ||.||i denotes the trace-norm flfill . In other words, in the 
limit of very small two-mode squeezing the maximally entan- 
gled state can be obtained to a high degree of accuracy. So 
the appropriate choice for the beam splitters on one side does 
depend on the value of q, whereas the beam splitter on the 
other side becomes redundant. In a similar manner, one can 
generate states of the form |0, 0) + axi|l, 1)- If one does not 
care about the phase of a\ t \, the correct choice for the above 
transmittivities and reflectivities is then 



l*(«)l 



«13i 



a 



(0)|2 
1,1 1 



1/2 



(59) 



and \r(q) \ := [1 — |£(g)| 2 ] 5. This analysis shows that with the 
help of passive optical elements and photon detectors, quan- 
tum states of the appropriate kind can in fact be prepared. 
There is, however, a trade-off concerning accuracy of the pro- 
tocol and success probability: For any finite q, the resulting 
states are not exactly pure, whereas the probability of success 
(such that the non-vacuum outcome is obtained in both detec- 
tors) is a monotone decreasing function of q. 

The resulting states of this protocol can then form the start- 
ing point of the generation of Gaussian states via the protocol 
in Sec. [H] In effect, this scheme allows one to generate ap- 
proximate Gaussian states (in fact, two-mode squeezed vacua) 
with a higher q than the initial supply, which is nothing other 
than a distillation procedure. 

An example of the results of such a distillation protocol, 
where the initial step is followed by three iterations of the pro- 
tocol from Sec. ||1J is illustrated in Fig. 0] The overall prob- 
ability is far lower than for three steps of the protocol from 
Sec.|H]alone (cf. Fig.|2j, due to the low success probability 
of the initial step. This is largely due to the low probability 
of measuring the presence of photons on the side where no 
beam splitter is employed, i.e. Alice's side. Since the effect 
of this measurement is to prepare a single photon on Bob's, 
this low probability step could be avoided if a single photon 
source were available. 

In light of the fact that distillation with Gaussian opera- 
tions alone was shown to be impossible 11(1 ll 111 , it is then 
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FIG. 4: This figure illustrates a full distillation procedure. Begin- 
ning with a supply of two-mode squeezed vacua, with q = 0.01, the 
protocol outlined in Sec. UVl is then applied, which maps this state 
onto a non-Gaussian state of higher entanglement, followed by three 
iterations of the protocol described in Sec. [2] The properties of state 
produced depend on the transmittivity T of the beam splitter em- 
ployed in the first step. Here, the factor by which the entanglement 
of the final achieved state -Eanai is greater than the entanglement of the 
initial supply E[ nt , (where the entanglement is calculated as the Von 
Neumann entropy of the reduced density matrix of a single mode) 
is plotted as a solid line, and the overall success probability of the 
entire process, when this initial step is followed by three iterations of 
the protocol to generate Gaussian states, is plotted as a dashed line. 



significant that this scheme does, in fact, realise pure-state 
distillation into approximate Gaussian states via suitable non- 
Gaussian operations, here photon detection. 

This simple protocol is not suitable when the initial sup- 
ply consists of two-mode squeezed states with a high q, and 
another method of generating non-Gaussian states of higher 
entanglement must be used. A more detailed analysis of op- 
timal preparation protocols that only include passive optical 
elements and photon detectors will be investigated elsewhere. 
Here, we concentrate on the proof-of-principle that Gaus- 
sian states can indeed be distilled to approximately Gaussian 
states. 



V. DISCUSSION AND CONCLUSIONS 

We have shown that, using passive optical elements and 
photon detectors that do not distinguish different photon num- 
bers, one can distill pure Gaussian states to arbitrarily high 
precision, in spite of the impossibility of distilling Gaussian 
states with Gaussian operations 11(1 ll ill . It should be noted 
that in our discussion we have assumed the photon detectors to 
have unit efficiency, in order to show that how one can in prin- 
ciple generate Gaussian states from a non-Gaussian supply. 
Needless to say, in any experimental realisation, one would 
have to deal with detector efficiencies significantly less than 
one. Such detectors can, e.g., be modeled by employing per- 
fect detectors, together with an appropriate beam splitter with 
an empty input port llisil . If the detector efficiency is still close 
to one, one would expect - after a small number of iterations 
of the procedure - the resulting states to be still close to those 
presented in this idealised protocol. The convergence proper- 
ties will in general be different from the ideal situation. Dark 
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counts of the detector, in turn, do not affect the performance of 
the protocol, except that the success probability is decreased. 
These matters will be discussed in more detail elsewhere. 

In several practical applications of the procedure, one can 
actually assume the initial state to be known. This is the case, 
for example, if one uses the above protocol in order to pu- 
rify a state in a quantum privacy amplification procedure 1 9 ] . 
Then, one may use homodyne detection together with passive 
optical elements in order to implement a positive-operator- 
valued measure (POVM) {|a)(a| : a £ (D}, where \a) de- 
notes the state vector of a coherent state, instead of photon 
detection 1 1 1, 20]. This would render a displacement in phase 
space necessary in the last step, depending on the measure- 
ment outcomes in each step. Such a modification, however, 
would transform the originally probabilistic protocol into a 
deterministic one. Also, the detector efficiencies can be as- 
sumed to be significantly larger. Even the displacement could 
be accounted for in the classical analysis of the measured data 
in the final stage of a protocol that makes use of the prepared 
entangled Gaussian state, e.g., a quantum cryptography proto- 
col. 

In this article we have restricted our analysis to pure states. 
In practical implementations it would clearly also be useful 
to be able to distill highly entangled Gaussian staes from a 
mixed initial supply. However, the full treatment of these pro- 



tocols for general mixed states is lengthy and will be presented 
elsewhere. To summarize, we have identified a procedure, 
that asymptotically produces Gaussian states from a supply 
of non-Gaussian, finite-dimensional states by means of Gaus- 
sian operations . In fact, the limiting Gaussian state for a pure 
given input can be found analytically. We have seen that even 
after a very small number of iteration steps the degree of over- 
lap between the resulting state and the theoretical limit state 
is close to unity. Moreover, the probability of obtaining this 
approximate state is of the order of 0.1. In that respect the 
whole protocol is experimentally feasible with present tech- 
nology. This result should contribute to the search for strate- 
gies to distribute continuous-variable entanglement over large 
distances. 



Acknowledgments 

We would like to thank I. Walmsley and J.I. Cirac for fruit- 
ful discussions, and K. Audenaert, A. Lund and W.J. Munro 
for very helpful remarks on the manuscript. This work was 
partially funded by the Alexander von Humboldt foundation, 
Hewlett-Packard Ltd., the EPSRC and the European Commis- 
sion (EQUIP). 



[1] R.E. Slusher, L.W. Hollberg, B. Yurke, J.C. Mertz, and J.F. Val- 
ley, Phys. Rev. Lett. 55, 2409 (1985); L.A. Wu, H.J. Kim- 
ble, J.L. Hall, and H. Wu, Phys. Rev. Lett. 57, 2520 (1986); 
M.D. Reid and D.F. Walls, Phys. Rev. A 34, 1260 (1986); 
Ch. Silberhorn, P.K. Lam, O. Weiss, F. Konig, N. Korolkova, 
and G. Leuchs, Phys. Rev. Lett. 86, 4267 (2001); C. Schori, 
J.L. S0rensen, and E.S. Polzik, Phys. Rev. A 66, 033802 (2002); 
M.M. Wolf, J. Eisert, and M.B. Plenio, Phys. Rev. Lett. 90, 
047904 (2003). 

[2] Gaussian states of systems with canonical degrees of freedom 
are those states with a Gaussian Wigner function (or equiva- 
lently, a Gaussian characteristic function) 11511 . 

[3] L. Vaidman, Phys. Rev. A 49, 1473 (1994); S.L. Braunstein and 
H.J. Kimble, Phys. Rev. Lett. 80, 869 (1998); G.J. Milburn and 
S.L. Braunstein, Phys. Rev. A 60, 937 (1999). 

[4] S.L. Braunstein and H.J. Kimble, Phys. Rev. A 61, 042302 
(2000). 

[5] T.C. Ralph and P.K. Lam, Phys. Rev. Lett. 81, 5668 (1998); 
M. Hillery, Phys. Rev. A 61, 022309 (2000); D. Gottesman 
and J. Preskill, Phys. Rev. A 63, 22309 (2001); F. Grosshans 
and P. Grangier, Phys. Rev. Lett. 88, 057902 (2002); P. Navez, 
Eur. Phys. J. 18, 219 (2002); Ch. Silberhorn, T.C. Ralph, 
N. Lutkenhaus, and G. Leuchs, Phys. Rev. Lett. 89, 167901 
(2002). 

[6] S.J. van Enk, J.I. Cirac, and P. Zoller, Science 279, 205 (1998); 

S. Scheel and D.-G. Welsch, Phys. Rev. A 64, 063811 (2001); 

S. Scheel and D.-G. Welsch, quant-pli/0207114. 
[7] CH. Bennett, D.P DiVincenzo, J.A. Smolin, and W.K. Woot- 

ters, Phys. Rev. A 54, 3824 (1996); see also Refs. Q|. 
[8] J.-W. Pan, C. Simon, C. Brukner, and A. Zeilinger, Nature 410, 

1067 (2001); X.B. Wang, quant-ph/0208 lES] 
[9] H. Aschauer and H.J. Briegel, Phys. Rev. Lett. 88, 047902 



(2002); H. Aschauer and H.J. Briegel, Phys. Rev. A 66, 032302 
(2002). 

[10] J. Eisert, S. Scheel, and M.B. Plenio, Phys. Rev. Lett. 89, 
137903 (2002). 

[11] J. Fiurasek, Phys. Rev. Lett. 89, 137904 (2002); G. Giedke and 
J.I. Cirac, Phys. Rev. A 66, 032316 (2002). 

[12] B. Demoen, P. Vanheuverzwijn, and A. Verbeure, Lett. Math. 
Phys. 2, 161 (1977); B. Demoen, P. Vanheuverzwijn, and 
A. Verbeure, Rep. Math. Phys. 15, 27 (1979); J. Eisert and 
M.B. Plenio, Phys. Rev. Lett. 89, 097901 (2002); J. Fiurasek, 
Phys. Rev. A 66, 012304 (2002). 

[13] L.-M. Duan, G. Giedke, J.I. Cirac, and P. Zoller, Phys. 
Rev. Lett. 84, 2722 (2000); T. Opatrny, G. Kurizki, and D.- 
G. Welsch, Phys. Rev. A 61, 032302 (2000); L.-M. Duan, 
G. Giedke, J.I. Cirac, and P. Zoller, Phys. Rev. A 62, 032304 
(2000); G. Giedke, L.-M. Duan, J.I. Cirac, and P. Zoller, Quant. 
Inf. Comp. 1(3), 79 (2001); PT. Cochrane, T.C. Ralph, and G.J. 
Milburn, Phys. Rev. A 65, 062306 (2002). 

[14] Note that this convergence is pointwise convergence for the co- 
efficients specifying the state vector. 

[15] W. Vogel, D.-G. Welsch, and S. Wallentowitz, Quantum Optics, 
An Introduction (Wiley-VCH, Weinheim, 2001). 

[16] R.A. Horn and C.R. Johnson, Matrix Analyis (Cambridge Uni- 
versity Press, Cambridge, 1987). 

[17] X. Ma and W. Rhodes, Phys. Rev. A 41, 4625 (1990). 

[18] S.D. Bartlett and B.C. Sanders, Phys. Rev A 65, 042304 (2002). 

[19] K. Nemoto and S.L. Braunstein, Phys. Rev. A 66, 032306 
(2002). 

[20] U. Leonhardt, Measuring the quantum state of light (Cambridge 
University Press, Cambridge, 1997). 



